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Abstract 

In this paper we develop a tractable structural model with analytical default 
probabilities depending on some dynamics parameters, and we show how to cal- 
ibrate the model using a chosen number of Credit Default Swap (CDS) market 
quotes. We essentially show how to use structural models with a calibration 
capability that is typical of the much more tractable credit-spread based in- 
tensity models. We apply the structural model to a concrete calibration case 
and observe what happens to the calibrated dynamics when the CDS-implied 
credit quality deteriorates as the firm approaches default. Finally we provide a 
typical example of a case where the calibrated structural model can be used for 
credit pricing in a much more convenient way than a calibrated reduced form 
model: The pricing of counterparty risk in an equity swap. 
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1 Introduction 

Modelling firms default is an important issue, especially in recent times where the 
market is experiencing an increasing interest in credit derivatives trading. Pricing 
models can be divided into two main categories: (i) reduced form models and (ii) 
structural models. 

Reduced form models (also called intensity models when a suitable context is 
possible) describe default by means of an exogenous jump process; more precisely, 
the default time is the first jump time of a Poisson process with deterministic or 
stochastic (Cox process) intensity. Here default is not triggered by basic market 
observables but has an exogenous component that is independent of all the default 
free market information. Monitoring the default free market does not give complete 
information on the default process, and there is no economic rationale behind default. 
This family of models is particularly suited to model credit spreads and in its basic 
formulation is easy to calibrate to Credit Default Swap (CDS) data. See references 
at the beginning of Chapter 8 in Bielecki and Rutkowski (2001) for a summary of 
the literature on intensity models. We cite here Duffie and Singleton (1999) and 
Lando (1998), and Brigo and Alfonsi (2003) as a reference for explicit calibration of 
a tractable stochastic intensity model to CDS data. 

Structural models are based on the work by Merton (1974), in which a firm life is 
linked to its ability to pay back its debt. Let us suppose that a firm issues a bond 
to finance its activities and also that this bond has maturity T. At final time T, 
if the firm is not able to reimburse all the bondholders we can say that there has 
been a default event. In this context default may occur only at final time T and 
is triggered by the value of the firm being below the debt level. In a more realistic 
and sophisticated structural model (Black and Cox (BC) (1976), part of the family 
of first passage time models) default can happen also before maturity T. In first 
passage time models the default time is the first instant where the firm value hits 
from above either a deterministic (possibly time varying) or a stochastic barrier, 
ideally associated with safety covenants forcing the firm to early bankruptcy in case 
of important credit deterioration. In this sense the firm value is seen as a generic asset 
and these models use the same mathematics of barrier options pricing models. For a 
summary of the literature on structural models, possibly with stochastic interest rates 
and default barriers, we refer for example to Chapter 3 of Bielecki and Rutkowski 
(2001). It is important to notice that structural models make some implicit but 
important assumptions: They assume that the firm value follows a random process 
similar to the one used to describe generic stocks in equity markets, and that it 
is possible to observe this value at any time. Therefore, unlike intensity models, 
here the default process can be completely monitored based on default free market 
information and comes less as a surprise. However, structural models in their basic 
formulations and with standard barriers (Merton, BC) have few parameters in their 
dynamics and cannot be calibrated exactly to structured data such as CDS quotes 
along different maturities. 

In this paper we plan to effectively use a structural model in the ideal "territory" 
of intensity models, i.e. to describe default probabilities in a way that is rich enough 
to calibrate CDS quotes. First of all, in Section [2] we show how standard barrier op- 
tions formulas can be used to find these probabilities, and in particular we investigate 
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the case of time dependent parameters in the dynamics of the underlying process. In 
Section [3] we show how to cahbrate the parameters of the model to market prices of 
single name CDS's (in a way similar to the procedure used to calibrate intensities in 
reduced form models). In Section H] we consider the concrete case of a firm approach- 
ing default, and we show how our proposed structural model calibration changes as 
the company credit quality (as summarized by its CDS quotes) deteriorates in time. 
The calibrated model is useful for subsequent pricing of more sophisticated deriva- 
tives depending on default. Indeed, in Section [5] we see an important example of how 
our structural model can be used to price some products involving considerations 
from both the credit and the equity market. The product we consider is an equity 
return swap with counterparty risk. In this context we have to take care of the cor- 
relation between the counterparty and the underlying, and this is done much more 
conveniently in structural models than in intensity models. 

Finally, we have considered tractable extensions of the ATIP model based on 
introducing scenarios on the value of the firm volatility and safety barrier. This 
framework allows to consider situations where the balance sheets information can be 
uncertain or hid some developments. This work and the related SVBATIP models 
are reported in Brigo and Tarenghi (2005), with a calibration case study on CDS 
data. 

2 The Black Cox approach and Barrier options 
formulas 

The fundamental hypothesis of the model we resume here is that the underlying 
process is a Geometric Brownian Motion (GBM), which is also the kind of process 
commonly used for equity stocks in the Black Scholes model. 

Classical structural models (Merton, Black Cox) postulate a GBM (Black and 
Scholes) lognormal dynamics for the value of the firm V. This lognormality assump- 
tion is considered to be acceptable. Crouhy et al (2000) report that "this assumption 
[lognormal V] is quite robust and, according to KMVs own empirical studies, actual 
data conform quite well to this hypothesis.". 

In these models the value of the firm V is the sum of the firm equity value S and 
of the firm debt value D. The firm equity value S, in particular, can be seen as a 
kind of (vanilla or barrier-like) option on the value of the firm V. Merton typically 
assumes a zero-coupon debt at a terminal maturity T. Black Cox assume, besides a 
possible zero coupon debt, safety covenants forcing the firm to declare bankruptcy 
and pay back its debt with what is left as soon as the value of the firm itself goes 
below a "safety level" barrier. This is what introduces the barrier option technology 
in structural models for default. 

More in detail, in Merton's model there is a debt maturity T, a debt face value L 
and the company defaults at final maturity (and only then) if the value of the firm 
Vf is below the debt L to be paid. 

The debt value at time t < T is thus 

A = MD{t,f) mm{Vf,L)] = Et[D{t, f)[Vf - {Vf - L)+]] = 
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= Et[D{t,T)[L - (L - Vf)+]] = P{t,T)L - Put(t,f ; Vt, L) 

where Put (time, maturity, underlying, strike) is a put option price, and the stochastic 
discount factor at time t for maturity T is denoted by D{t,T) = B{t)/B{T), where 
B(t) = exp^J^ Tudu) denotes the bank-account numeraire, r being the instantaneous 
short interest rate. Since we will assume deterministic interest rates, in our case 
D{t,T) = P{t,T), the zero coupon bond price at time t for maturity T. 

The equity value can be derived as a difference between the value of the firm and 
the debt: 

St = Vt-Dt = Vt- Pit, f)L + Put(t, f ; Vt, L) = Call(t, f ; Vt, L) 

so that, as is well known, in Merton's model the equity can be interpreted as a 
call option on the value of the firm. 

Let us now move to the Black Cox (BC) model. In this model we have safety 
covenants in place, in that the firm is forced to reimburse its debt as soon as its value 
Vt hits a low enough "safety level" H{t). The choice of this safety level is not easy. 
Assuming a debt face value of L at final maturity T as before, an obvious candidate 
for this "safety level" is the final debt present value discounted back at time t, i.e. 
LP{t,T). However, one may want to cut some slack to the counterparty, giving it 
some time to recover even if the level goes below LP{t, T), and the "safety level" can 
be chosen to be lower than LP{t,T). 

We will come back to this issue later on. For the time being we just assume V to 
be a geometric Brownian motion. 

dv; = v; jj,; dt + v; at dw{t), Vq = Vo (i) 

and impose a time dependent safety barrier H*{t). For the time being let us assume 
to have constant parameters (i.e. fj,^ = fj,* and at = a) and a constant barrier 
{H*{t) = H*). We set ourselves in the risk neutral measure, so that the drift term 
has the standard form of interest rate minus a payout ratio (i.e. /i* = r — q*) and 
([1]) is a standard Geometric Brownian Motion. Let us pretend for a moment that 
we are just doing derivatives pricing in an equity market with underlying V. In the 
literature it is possible to find many analytical formulas for barrier options pricing 
such as, for example, knock-in or knock-out options, and also digital options. One 
particular case is the down and out digital option or down and out bond (DOB), 
that is a contract paying one unit of currency at maturity T if, between the starting 
date of the contract and its maturity T, the underlying never touches the barrier H* 
{H* < Vq) from above. If we go back to the default interpretation with H* as safety 
barrier level and T < T, this option price is actually the price of a defaultable zero 
coupon bond with no recovery in the structural model framework. Now, if we call r 
the first time instant where the process hits the barrier from above, the price of this 
option in the risk-neutral framework is given by 

DOB(0,T)=E{D(0,T)l{.>r}} (2) 

where E is the risk neutral expectation associated with the risk neutral measure Q. 
Under our deterministic rates framework (or under interest rates that are independent 
of the underlying process) we can write 

DOB(0, T) = P(0, T) E{1|.>T}} = P(0, T) Q{r > T}. (3) 
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The last factor is the (risk-neutral) probability of never touching the barrier before 
T, also called survival probability. By means of stochastic calculus it is possible to 
derive analytically the price of the option by explicitly computing ([3]) (see for example 
Bielecki and Rutkowski (2001)): 



DOB(0,T) = P(0,T) 



$ (rfi) - 



H* 



l-2(r-g')/(<x2) 



$(4) 



(4) 



with di^2 = (^ilog;^ + (^r — q* — /(^"V^). By comparing equations ([3]) and 



we obtain 



It>T} 



l-2(r-q*)/(a2) 



$(c/2 



(5) 



Similar formulas can be obtained also for a barrier that is not necessarily constant 
in time (or "flat"), but with a particular exponential shape (see again Bielecki and 
Rutkowski (2001) for more details). 

When relaxing the assumption of constant parameters in the V* dynamics, the 
situation becomes much more complicated. In this case, even with a fiat barrier, it is 
not possible to find closed form pricing formulas. However, some recent work dealing 
with option pricing on underlying assets having time dependent parameters in the 
dynamics shows that it is possible to find analytical barrier option prices when the 
barrier has a particular curved shape depending partly on the dynamics parameters. 
See for example Rapisarda (2003), who builds on the fundamental work of Lo et al. 
(2003), where the formulas are expressed in a shape resembling the classical constant 
coefficients Black-Scholes formulas. Again, let us assume an underlying process like 
([1]), with = Tt — ql, where r and q* are respectively the time varying instantaneous 
risk free interest rate and payout ratio, and consider an option with maturity T. Let 
us take a barrier of the form: 



H*{t) =Hexp 



-q:-{l + 2f3) 



ds 



(6) 



depending on a parameter (3 and on the constant reference value H. For this special 
barrier it is possible to obtain exact formulas for many different barrier options when 
the underlying is V*. Obviously, if we had to price an option with a given constant 
barrier, the price obtained with these formulas would not be properly correct, but 
it would be a good approximation if we chose the barrier profile H*{t) as close as 
possible to the given fiat barrier value. It can be shown that the value of (3 keeping 
the barrier as fiat as possible is given by 



(3* 



and this reduces to P* = {r — q* 
(where, as expected, H*{t) = H). 



^ I ds 



if' 



aids I dt 



aMs ] dt 



(7) 



(T^/2)/(cr^) in the case of constant coefficients 
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Let us go back to the time- varying barrier 1^. In this framework 



DOB(0,T) = P(0,T)- 



Vo J 



where Vt = rt — ql — ^ and, using again ([2]), we find 



[r>T} 



log t + /o^ -^ds \ _ /i/- (0) y% ( log ^ + vsds 



lo ^sds 



Vo 



So "^sds 



where we recall that H*{0) = H exp (^r, - q* - (1 + 2/3)4) ^■^) • 

In the next section we show how this formula can be used to calibrate a given 
structural model to market data, in particular to CDS's quotes. 



3 Calibration of the structural model to CDS data 

Since we are dealing with default probabilities of firms, it is straightforward to think of 
financial instruments depending on these probabilities and whose final aim is to pro- 
tect against the default event. One of the most representative protection instruments 
is the Credit Default Swap (CDS). CDS's are contracts that have been designed to 
offer protection against default. Consider two companies "A" (the protection buyer) 
and "B" (the protection seller) who agree on the following. 

If a third reference company "C" (the reference credit) defaults at a time tq G 
{Ta,Tb], "B" pays to "A" at time r = tq itself a certain "protection" cash amount 
Lgd (Loss Given the Default of "C" ), supposed to be deterministic in the present 
paper. This cash amount is a protection for "A" in case "C" defaults. A typical 
stylized case occurs when "A" has bought a corporate bond issued from "C" and is 
waiting for the coupons and final notional payment from this bond: If "C" defaults 
before the corporate bond maturity, "A" does not receive such payments. "A" then 
goes to "B" and buys some protection against this risk, asking "B" a payment that 
roughly amounts to the bond notional in case "C" defaults. 

Typically Lgd is equal to a notional amount, or to a notional amount minus a 
recovery rate. We denote the recovery rate by "Reg" . 

In exchange for this protection, company "A" agrees to pay periodically to "B" 
a fixed "running" amount R, at a set of times {Ta+i, . . . , Tf,}, = Tj — Tj_i, Tq = 0. 
These payments constitute the "premium leg" of the CDS (as opposed to the Lgd 
payment, which is termed the "protection leg"), and R is fixed in advance at time 0; 
the premium payments go on up to default time r if this occurs before maturity Tf,, 
or until maturity Tf, if no default occurs. 

"B" protection Lgd at default tc if < tc < U "A" 

"B" ^ rate R at T„+i, . . . , T,, or until default tc ^ "A" 

Formally, we may write the RCDS ("R" stands for running) discounted value at 
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time t seen from "A" as 

b 

UliCVSa,bit) := ~D{t,T){T -Ti3(^r)~l)R'i-{Ta<T<Ti} - ^ D{t,Ti)aiRl{r>T,} 

i=a+l 

+ MT^<r<n}D{t,T) Lgd (10) 

where t G [T/3(j)_i, T^(t)), i.e. Tpi^t) is the first date among the Tj's that follows t, and 
where is the year fraction between Tj_i and Tj. 

Sometimes a slightly different payoff is considered for RCDS contracts. Instead 
of considering the exact default time r, the protection payment Lgd is postponed to 
the first time Tj following default, i.e. to Tg(r)- If the grid is three or six months 
spaced, this postponement consists in a few months at worst. With this formulation, 
the CDS discounted payoff can be written as 

b b 

+ l{T^^^<r<T,}D{t,Ti) Lcn, (11) 

i=a+l i=a+l 

which we term "Postponed Running CDS" (PRCDS) discounted payoff. Compare 
with the earlier discounted payout ffTOl) where the protection payment occurs exactly 
at r: The advantage of the postponed protection payment is that no accrued-interest 
term in (r — T^(t-)_i) is necessary, and also that all payments occur at the canonical 
grid of the Tj's. The postponed payout is better for deriving market models of CDS 
rates dynamics and for relating CDS's to floaters, see for example Brigo (2004, 2004b). 
When we write simply "CDS" we refer to the RCDS case. 

Let us consider again the basic RCDS: The pricing formula for this payoff depends 
on the assumptions on the interest rates dynamics and on the default time r. Let 
J-'t denote the basic flltration without default, typically representing the information 
flow of interest rates and possibly other default-free market quantities (and also in- 
tensities in the case of reduced form models), and = J^t V a ({r < u},u < t) the 
extended flltration including explicit default information. In our current "structural 
model" framework with deterministic default barrier the two sigma-algebras coin- 
cide by construction, i.e. Qt = ^t, because here the default is completely driven by 
default-free market information. This is not the case with intensity models, where 
the default is governed by an external random variable and J^t is strictly included in 
Gt, i-e. J^t C Gf 

We denote by CDS(t, [Ta+i, . . . , Ti,],Ta, Ti,, R, Lgd) the price at time t of the above 
standard running CDS. At times some terms are omitted, such as for example the list 
of payment dates [T^+i, . . . , Tf,]. In general we can compute the CDS price according 
to risk-neutral valuation (see for example Bielecki and Rutkowski (2001)): 

CDS{t,Ta,Th,R, Lgd) = E{UKCUSa,b{t)\gt} = E{UKCUSa,b{t)\J^t} =: Et{UKCUSa,b{t)} 

(12) 

in our structural model setup. A CDS is quoted through its "fair" R, in that the rate 
R that is quoted by the market at time t satisfles CDS(t, T^, T;,, i?, Lgd) = 0. Let us 
assume, for simplicity, deterministic interest rates; then we have 

CD5(t,T,,T,,i?, Lgd) := -i?Ei{P(t,r)(r-T^(.)_i)l|T„<r<T,}} 

b 

{Ta<T<n}Pit,T)}. (13) 

i=a+l 
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It is clear that the fair rate R strongly depends on the default probabilities. The idea 
is to use quoted values of these fair i?'s with different maturities to derive the default 
probabilities assessed by the market. 

While in simple intensity models the survival probabilities can be interpreted as 
discount factors (with credit spreads as discounting rates), and as such can be easily 
stripped from CDS's or corporate bonds, in structural models the situation is much 
more complicated. In fact, here, it is not possible to find a simple "credit spread" 
formulation for dQ{T > t} starting from (Q. Moreover we have to pay attention 
to one fundamental aspect of the model. The barrier triggering default depends on 
the maturity chosen for a particular instrument, i.e. H*{t) = H^{t) has a parametric 
dependence on T. But then, if we plan to use the model for instruments with different 
maturities, it is necessary to impose a consistency condition, that is 

H^,it) = H^t) (14) 

for t < Ti < T2 for every pair of maturities Ti and T2. With our choice of H*{t) 
above we have easily that this condition implies, when enforced for all possible Ti, T2: 

2 2 
rt - ql - (1 + 2/3)^ = 0, or g* = r, - (1 + 2(3)^. (15) 

We notice that if this condition is satisfied, we have that the barrier if*, defined 
in terms of g*, obviously flattens to the constant value H. Thus consistency of 
the curved harrier for all maturities induces a flat harrier. In general, however, we 
expect (|T5|) not to hold if V* is the value of the firm, since r, q* and a are given 
to us exogenously and (3 is just constant. However, we can still manage to preserve 
analytical tractability as follows. Let us assume from now on that the real risk neutral 
dynamics of the firm value is a process given by 

dVt = Vt {rt - qt) dt + Vt dW{t) (16) 

where qt is the true payout ratio and q^ is defined by f fTSj) . The "system" (y*,H*) 
with said q* is both tractable and satisfying the consistency condition. Our problem 
is that V* is not the real firm value, so that the system {V*,H*) is not good for 
modeling default. However, if we define 

H{t) = Hexp j\q, - q^ds^ (17) 

by integrating V* and V^s equations it is easy to show that the first time V^* hits 
H*{t) = H is the same as the first time the real process Vt hits H{t). Therefore, 
default probabilities computed with the tractable model V*,H* are the same as 
default probabilities for the "true" model V, H. We can compute quantities and 
perform our calibration with the former model and consider the latter as the real 
model. We can thus state the following 

Proposition 3.1. (Analytically- Tractable First Passage (ATIP) Model) As- 
sume the risk neutral dynamics for the value of the firm V is characterized hy a 
risk free rate rt, a payout ratio qt and an instantaneous volatility at, according to 
equation i-G. 

dVt = Vt{rt - qt) dt + VtatdW{t) 



D. Brigo, M. Tarenghi: CDS Calibrated Structural Model and Application to Equity Swaps 10 



and assume a default barrier H{t) of the form given in [T^ with q* given as in l[T5\) . 
i.e. ^ ^ 

Hit) = Hexp ^ (^q, - r, + (1 + 2/3)^^ ds^ 

and let t be defined as the first time where V hits H from above, starting from Vq > H, 

r = mf{t >0:Vt< H{t)}. 
Then the survival probability is given analytically by 



> T} 



/o ^sds I \ \J lo ('Ids 



Formula ( fTSi) is easily obtained by substituting the consistency condition ( |T5i) in 
formula ([9]) for the survival probability. From our earlier definitions, straightforward 
computations lead to the price at initial time of a CDS, under deterministic interest 
rates, as 

cmaM Lgd) =R r ^(0, t){t - Tp^t)^^)dq{T > t) (19) 

~R J] P(O,T0aiQ(r > T,) - Lgd j " P{0,t)dQ{T > t) 

i=a+l '''^'^ 

SO that if one has a formula for the curve of survival probabilities t i— >■ Q(r > t), as 
in our ATIP structural model, one also has a formula for CDS. 

Notice an important feature of Formula (fT8|) and then (fT9!) : Survival (and default) 
probabilities and CDS values only depend on the ratio between V and H, and not on 
V and H separately. This means that, as far as default probabilities are concerned, a 
precise estimation of V is not needed. We can express default probabilities in relative 
terms, i.e. in terms of V/H. In other terms, we can imagine to re-scale V (considering 
V/Vq so that the initial condition reads Vq/Vo = 1) and also the barrier parameter H 
(taking H/Vq, a number smaller than one as barrier). 

Formula f|T8l) can be used to fit the model parameters to market data. However, 
the only parameter left that can account for time dependence is the volatility. If we 
use exogenous volatility (deduced perhaps from historical or implied equity volatility) 
we are left with no freedom. However, we may infer the first year volatility a{0-7-ly) : = 
{cr{t) : t E [0, 1]} from equity data and use if as a first fitting parameter, and then use 
the remaining later volatilities a{ly-r2y), a{2y-^3y) etc as further fitting parameters. 
Therefore a{2y^-) will be determined by credit quality as implied by CDS data rather 
than by equity data. To sum up, we can choose piecewise constant volatility, and look 
for those volatility values after the first year that make the quoted CDS's fair when 
inserting in their premium legs the market quoted i?'s. In this way we find as many 
volatilities as many CDS's we consider minus one. In this first case H is determined 
by CDS's and we explain below how exactly we find this H (credit spread method). 

Alternatively, if we aim at creating a one to one correspondence to volatility 
parameters and CDS quotes, we can exogenously choose the value H (for example we 



D. Brigo, M. Tarenghi: CDS Calibrated Structural Model and Application to Equity Swaps 1 1 



will see below the protection/excursion method, but judgemental analysis or trial and 
error are often needed) and jS, leaving all the unknown information in the calibration 
of the volatility. If we do so, we find exactly one volatility parameter for each CDS 
maturity, including the first one. 

In general the above CDS calibration procedures are justified by the fact that 
in the end we are not interested in estimating the real process of the firm value 
underlying the contract, but only in reproducing risk neutral default probabilities 
with a model that makes sense also economically. While it is important that the 
underlying processes have an economic interpretation, we are not interested in sharply 
estimating them or the capital structure of the firm, but rather we appreciate the 
structural model interpretation as a tool for assessing the realism of the outputs 
of calibrations, and as an instrument to check economic consequences and possible 
diagnostics. 

Finally, in case we still have a preferred terminal maturity T for the debt, and an 
indication of the final debt at maturity, which we call L, it makes sense to impose 
that our barrier H be always below the present value of the final debt, i.e. P{t,T)L. 
This condition amounts to assuming that we are cutting some slack to the firm by 
allowing it to go somehow below the debt present value before forcing it to declare 
bankruptcy. How much below can be decided by means of H and /?. In detail, our 
condition reads 

H exp - r, + (1 + 2/9)^^ < exp ^- ^ r,c/s j L. 

We can easily rewrite this condition as 



L > H exp 

A sufficient condition for this, in case the round brackets term is positive, is 

H < LP(0,f), 

i.e. the final debt initial present value has to larger than the barrier parameter 
H. This means that in all our use of the model with positive round brackets one may 
presume a final debt L at a preferred maturity (typically larger than ten years, since 
this is the largest CDS maturity) satisfying this condition. 

We are still in need to connect equity and firm value if we aim at deriving part 
of the value of the firm volatility from equity data. Since default is enforced as soon 
as V hits H at time r if before T, or is given at T if is below L, the debt value at 
time t would be 

A = MD{t, r)if(r)l{,<f }] + Et[D{t, f) min(L, )l{r>f }] 

which can be computed as a function of Vt with computations similar to f|T8l) (with 
time t replacing time 0). In turn, we would have the equity value a.s St = Vt — Dt = 
Vt — DtiVt) from which we could derive an approximation for Vs (unknown) volatility 
in terms of S"s (known, be it historical or implied) volatility through Ito's formula 
and some approximations. 
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However, we will not pursue this strategy but simply take the equity volatility 
itself as a proxy for the order of magnitude of the firm value volatility. As a matter 
of fact, our preferred approach will be to let CDS data select most (and in some cases 
all) values of the firm volatilities. 

3.1 First Numerical Example 

In this section we present some results of the calibration performed with the structural 
model. We consider CDS contracts having the telecoms sector Vodafone company 
as underlying with recovery rate Rec = 40% (Lgd = 0.6). In Table [1] we report 
the maturities of the contracts and the corresponding "mid" CDS rates -R^™(0) 
(quarterly paid) on the date of March 10th, 2004, in basis points {Ibp = 10~^). We 
take Ta = in all cases. 



Maturity Tf, 


Rate R^f{0) (bps) 


March 21st, 2005 


21.5 


March 20th, 2007 


33.0 


March 20th, 2009 


43.0 


March 21st, 2011 


49.0 


March 20th, 2014 


61.0 



Table 1: Maturities of Vodafone's CDS's with their corresponding rates on March 10, 2004. 

In Tables [2] and [3] we present the results of the calibration performed with the 
structural model and, as a comparison, of the calibration performed with a deter- 
ministic intensity (credit spread) model (using piecewise linear intensity). In this 
first example the parameters used for the structural model have been selected based 
on qualitative considerations, and are P = 0.5 and H/Vq = 0.5. The a's have been 
found by calibration to CDS quotes. Below we report both the values of the cal- 
ibrated parameters in the two models (volatilities and intensities) and the related 
survival probabilities. 



Maturity Tf, 


Volatility cr{Tb_i ^ %) 


Survival Prob. Q(r > Tb) 


March 10th, 2004 


24.343% 


100.000% 


March 21st, 2005 


24.343% 


99.625% 


March 20th, 2007 


12.664% 


98.315% 


March 20th, 2009 


12.766% 


96.352% 


March 21st, 2011 


12.659% 


94.204% 


March 20th, 2014 


15.271% 


89.645% 



Table 2: Calibrated piecewise constant volatilities nodes and subsequent survival probabil- 
ities with the structural model. In the first column there are the maturities of the contracts 
(the first entry is the date where the calibration is performed. All cr's have been obtained 
from CDS quotes, H and j3 being fixed exogenously) . 

As a further comparison, we plot the behavior of the survival probabilities T i— > 
Q{r > T} resulting from CDS calibration in both cases (Figure [1]). It is clear that the 
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Maturity 


Intensity 


Survival Prob. Q(r > Tf,) 


March lOtn, 2004 


0.357% 


100.000% 


iViaiCll ZibL, zuuo 


U.OO ( /o 


QQ f^97% 


March 20th, 2007 


0.952% 


98.316% 


March 20th, 2009 


1.033% 


96.355% 


March 21st, 2011 


1.189% 


94.206% 


March 20th, 2014 


2.104% 


89.604% 



Table 3: Intensity model. In the second column we present the nodes of the piecewise 
linear calibrated intensities and in the last column the corresponding survival probabilities. 
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Figure 1: Survival probabilities T i-^ Q{t > T} comparison between the CDS calibration 
with the structural model and the CDS calibration with the intensity model. The dotted 
line corresponds to survival probabilities computed with the intensity model, the continuous 
line to survival probabilities computed with the structural model. 

survival probability found with the two different models is nearly the same, which 
indirectly says that CDS are actually instruments that efficiently translate default 
probabilities into prices, since the default probabilities depend little on the chosen 
model, as long as this is consistent with the same CDS quotes. This is somehow 
obvious when keeping in mind (fT9l) . 

A difference is in the first period, i.e. during the life of the first CDS, where 
the probability computed with the structural model is slightly above the other one, 
given the less sudden nature of the default in structural models in general. This is 
illustrated in Figure [TJ This leads us to the following 

Remark 3.2. (Short term credit spreads). It is often said that structural models 
imply unrealistic short-term credit spreads. However, if "short-term" is meant as a 
realistic short maturity, our model does not suffer from this drawback, since it can 
calibrate any realistic 6m or hypothetical 3m CDS quote. If one is not happy with 
the default probability between and ly, it suffices to calibrate a shorter term CDS, 
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and the model can do this exactly. What is more, we will see below a case with very 
high default probability where only the structural model is able to calibrate the CDS 
quotes, the intensity model giving negative intensities as outputs. 

In Figure [2] we plot the term structure of the cahbrated volatihty (see also Table [2]). 
Sometimes the calibration yields values that are not plausible as volatility levels, but 
these volatilities are to be taken as default probabilities fitting parameters with an 
economic interpretation rather than directly as firm value volatilities. In any case, 
consistency with CDS quotes for relative value pricing is an important aspect we need 
to incorporate in any relative value pricing model. Also, embedding CDS information 
in a structural model framework allows us to check a posteriori the realism of the 
calibrated quotes and of the model, so that we may play with our degrees of freedom 
(for example in (3 and H) to obtain more realistic calibrations. 



30.00% 
25.00% 



20.00% 
15.00% 



10.00% ^ 1 

10-Mar-04 09-Apr-06 08-May-08 07-Jun-10 06-Jul-12 05-Aug-14 



Figure 2: Volatility term structure T ^ cr(T) obtained from the calibration. The corre- 
sponding data are presented in Table [21 

We can double check the validity of our results above by means of a Monte Carlo 
simulation. We can simulate numerically the process V* and check whether it crosses 
the default barrier H* or not: In both cases we can compute the (discounted) payoff 
of the CDS for that particular path. We repeat this procedure N times and find the 
expected value of the payoff as a sample mean from the simulated scenarios. If, for 
each maturity we use the rate i?o,6(0) given by the market, then we should obtain 
CDS's with zero values. 

First of all we derive the right dynamics for the underlying: We use the process ([1]) 
using the consistency condition f|T5l) . obtaining 

2 

dv; = v; (1 + 2p)^dt + v; at dw{t). (20) 

We know that in general simulations of processes with barriers require a very large 
number of paths and a very small time step. To avoid this problem we applied the 
Brownian Bridge method (see for example Metwally and Atiya (2002)), which allows 
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US to use a larger time step (actually we used an interval of about five days), thus 
reducing significantly the simulation time without losing precision on the reproduced 
default probabilities. 

In Table m we present the result of the Monte Carlo simulation. We used the same 
data considered at the beginning of the section, and the rates R given in Table [H We 
used = 250000 scenarios, finding that the values of the CDS's simulated according 
to the analytically calibrated dynamics are practically zero as should be (zero is inside 
the small standard error MC window). 



Maturity 


Simulated Payoff 


MC std. dev. 


March 21st, 2005 


0.2 


0.7 


March 20th, 2007 


0.7 


1.5 


March 20th, 2009 


-0.9 


2.1 


March 21st, 2011 


-0.8 


2.5 


March 20th, 2014 


-0.1 


3.1 



Table 4: Simulated payoffs of the previously considered CDS's in basis points. The errors of 
the Monte Carlo simulations are reported as well. The number of scenarios is = 250000. 



3.2 Finding /f: Credit spread method and protection/excursion 
analogy method 

Not all possible values of the parameters H and (3 are good for the calibration. The 
parameter f3 can be used to shape the safety covenants barrier H once q, a and r are 
given. And if we calibrate the a's, the initial choice of (3 will have repercussions on 
the calibration output patterns. In general for a given j3, if H is too high (or too 
low) the calibration may not succeed in finding all the volatilities. A way to choose a 
valid reference value H having also a more direct link with market information is the 
following. Suppose we know the equity volatility for one-year maturity (for example 
an implied volatility taken from quoted options) and suppose we set the structural 
model's a up to one year to a constant value given by said equity volatility. Then 
the probability of not defaulting before T (= 1 year) is given by (see ([11 



lT>T} 



logH + Pa'T 



(we have re-scaled V to one, i.e. Vq = 1, so that the found H will have to be multiplied 
by Vq) and depends on H. Being all other parameters known, ijj is monotonically 
decreasing in H. But if we adopt a deterministic intensity model for one moment, 
the market provides us with 

i-T 



>T} = exp \cDs{.s)ds^ 



where Xcds are the deterministic intensities that have been stripped from CDS data. 
To impose a reasonable barrier H we solve 



exp XcDs{.s)ds^ = i){H] p, T, a) 
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in H. Since ip is monotonic, this should be easily dealt with numerically. We may 
call this value the "credit spread based H" , since H comes from a default probability 
based on credit spreads stripped by CDS's. Using this H, we calibrate the structural 
model volatilities to the CDS quotes, including the first one. Since the intensity-based 
ly default probability extracted from CDS and the CDS quote itself contain similar 
information but are not exactly equivalent, we will find a a{0-rly) that is close to our 
initial one-year a selected from the equity market, but not exactly the same. This is 
due for example to the first CDS having quarterly payments on whose dates survival 
probabilities depend on the chosen model. Even if two models agree on the one-year 
survival probability, they are not necessarily agreeing on the say semiannual default 
probability. And if a payoff involves said default probability the two models may give 
slightly different values even if they agree on the one-year default probability. This 
is why for example in Table [7] we find a(0 -j- ly) = 5.012% when the equity volatility 
is 5%. 

A different possibility for choosing H is the following. We may compare the range 
of protection Lgd offered with the CDS to a rough covering of the excursion V — H 
characterizing default in the structural model. If we consider then LgdVq as the 
protection (on the notional Vq associated with the initial value of the firm) and set 
it equal to Vq — if we find 

Lgd = 1 - H/Vo, or H/Vo = Reg, 

so that we have an immediate candidate for H. We term this value the "excur- 
sion/protection analogy H" . However, if we resort to this value we may have numer- 
ical problems in stripping volatilities and moreover we cannot impose the one-year 
volatility from the equity market any longer, since we have already used the degree of 
freedom provided by H. Indeed the calibration proceeds in the same way as before, 
but problems may arise in low recovery cases. Low values for recovery rates imply 
low barriers, which in turn lead to high values for stripped volatilities, even much 
larger than 100%. Under these conditions there are times when the solver is not able 
to retrieve a feasible solution for the calibrated volatilities. A way to bypass this 
obstacle is given by the parameter (3: its main role is to vary the steepness of the 
safety covenant barrier, as is clear from (1171) . and this fact has a relevant impact also 
on the term structure of volatilities resulting from the calibration, as we explain more 
in detail in the following section. As a matter of fact, it is possible to indirectly act 
on the volatilities by varying (3, looking for those values of (3 returning reasonable 
values of a from the calibration, while keeping an eye on the realism of the resulting 
safety barrier. 

4 Further Model Specification and a Case Study: 
the Parmalat Crisis 

We now apply the model to a concrete case. We are going to see what happens, 
in terms of model parameters, when a firm approaches default. In particular we 
consider the case of Parmalat Finanziaria SpA, an important Italian company whose 
main activity is the production of milk and other alimentary goods, and also present 
in a large variety of other sectors. In the past years this company and the related 
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stock have been very important for the Itahan economy; in fact in 2003 the Parmalat 
stock was included in the main Itahan stock index (Mib 30). After some signal 
worrying the investors, at the end of 2003 the Italian control commission (Consob) 
began to analyze deeply the balance sheets of the company, and as the investigation 
proceeded the company entered a deep crisis. 

Here we plan to analyze how the structural model behaves in such a situation, 
as the credit quality implied by CDS quotes deteriorates in time. First of all, in 
Figure [3] we present the graph of the stock price and the graph of the (ly-) CDS rate 
R. It is immediate to see that there exists a sort of negative correlation between the 
two, as expected: As the equity value S increases, the CDS premium R (expressing 
a protection cost) decreases. From the historical data, we see that the stock price 
experienced a deep fall on December 11th. So we consider some particular dates 
before that day, some in "standard regime conditions" and also some dates in "crisis 
conditions" , to see how the CDS calibration changes with market conditions. 



STOCK PRICE AND CDS SPREAD (ly) 
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Figure 3: Historical behavior of the Parmalat stock St (in euros) and of the ly-CDS R{t) 
(in bps). The dotted Une corresponds to the stock price (left scale) while the continuous 
Hne to the CDS rate (right scale). 

In Table |5] we report the considered dates, with the related stock prices, CDS 
rates (in basis points), and the recovery rates. We also add the volatilities: These 
are historical volatilities obtained by the last month of data (as records of implied 
volatility were not available). These volatilities have been sometimes approximated 
in excess, to take into account the fact that in crisis periods the stock becomes more 
volatile than in the previous period, an effect that is not modelled by historical 
volatility. 

Some important remarks on the data. First of all we see that the CDS rates R 
(which were already high in the first dates we considered) grow very fast and the 
recovery rate decreases as the crisis approaches. The historical volatility we use as a 
proxy for ly volatility increases a lot in time as well (probably the last value of 50% 
is even underestimated, considering that after a couple of days the price has gone 
below one euro). We have also set the payout ratio identically equal to zero. 
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14 
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Table 5: Input data. 

With this data, we are now able to compute the credit spread based barrier 
parameter H for each considered date. In principle we could plot the value H on 
a graph with the corresponding initial value Vq (either normalized or not), and also 
the behavior of the real barrier H. Even if we did this for each date, this would not 
be sufficient to deduce whether the model efficiently incorporates default information 
because the different volatilities may be misleading. More precisely, for normalized 
values about 1 we find H = 0.8977 for September 10th and H = 0.7253 for December 
10th. As the first safety barrier is closer to the initial value Vq = 1, one could think 
that the default is more likely in the first situation, but we know that this is not the 
case. The fact is that in the first situation we have a much lower volatility than in 
the second. So we decided to plot not only the barrier H and the initial value of the 
firm Vo, but also the expected price in the future t \—>- K{Vt) and a confidence region 
t ^ exp{E[\n{Vt)]-Std[\n{Vt)]) and t ^ exp(E[ln(l^t)] + Std[ln(l^i)]), where Std is the 
risk neutral standard deviation. Std is easily computed for normal variates such as 
In(l^). All these graphs are reported in Figures from H] to [71 In the confidence region 
between the expected price minus/plus one standard deviation there are the 68% of 
the total scenarios. Hence the instant in which the lower bound of the confidence 
region hits the barrier is a signal of the proximity of the default. And from the 
graphs it is clear that as time passes, the default becomes more and more likely and 
the scenarios trajectories widen, hitting the barrier earlier. Indeed, if on September 
10th the confidence region hits the barrier after nearly ten years, on December 10th 
this happens after less than one year. At this point of the discussion a few remarks 
on the choice of [3 are in order. 

Putting the consistency condition f[T^ in the equation for the barrier H given 
in (fTT]) . after some calculations we find 



H{t) = ^E[Vt]exp 







(21) 



Remark 4.1. The role of the (3 parameter. The integral of the squared volatility 
has the effect of introducing a concave and decreasing factor into the shape of the 
safety barrier, provided the round brackets term is positive. The role of (3 is thus clear: 
For given values of the volatilities a , (3 is the parameter that when increased allows 
us to increase the concavity and the steepness of the decreasing factor in the barrier, 
thus lowering it. This parameter may be useful when the calibrated volatilities get out 
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of control. Indeed, in case the needed concavity is high and (3 is small, volatilities may 
tend to explode to obtain the required harrier factor. By increasing [3 the volatilities 
needed to produce a given steepness and concavity can become smaller, avoiding this 
risk. 

In our first experiments we choose (3 = 0.5 such that the argument of the expo- 
nential is simply the integral of the variance, but this is just our particular choice 
based on numerical experiments. In this case, for low values of the volatilities the 
barrier is very similar to the behavior of the expected price, hence it is upward sloped. 
On the contrary, high values of the volatilities change significantly the shape of the 
barrier, but also the confidence region widens as a consequence. These effects can be 
checked by comparing Figures from H] to [7] to the Tables reported in the Appendix. In 
Figure [S] we show a typical term structure of volatility obtained from the calibration, 
corresponding to the data of December 8th. 



1.6 



1.4 



1.2 



1 




0.8 ^ < < < < 1 

10-Sep-03 09-Oct-05 08-Nov-07 07-Dec-09 06-Jan-12 04-Feb-14 



Figure 4: Default barrier t ^ H{t) on September 10th, 2003, t ^ E{Vt) and t 
exp (E[ln(I4)]^Std[ln(yf)]). 

Furthermore, in Figure [9] we plot the survival probability on December 10th. 
We see that the probability falls down very quickly in the first period, but later it 
decreases in a much slower way. This fact can be explained as follows. When the 
firm is approaching default, investors are more interested in buying protection over a 
small time horizon, because they assess a very high probability to the default event 
for the immediate future period. But if the firm survives to the full crisis, probably it 
will reach a certain stability situation, characterized by a "lower" default probability, 
so that provided the company survives the first immediate times after the crisis, 
investors are less worried by subsequent periods. 

As a final investigation, we have repeated the calibration when H is selected under 
the alternative method of the "excursion/protection analogy", i.e. by setting H/Vq = 
Rec. In the Appendix we report all the results. The first important remark is that here 
we have used (3 = 0.08. This is due to the fact that the calibration for the last date 
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Figure 6: Same as above on December 8th, 2003. 



requires very high volatihties, being the relevant H very low. Then we are compelled 
to use low /3's to avoid the calibration explosion hinted at in Remark 14.11 At the 
same time we cannot diminish f3 too much, because the second volatility o'{2y^3y) is 
very low due to the peculiar structure of CDS market quotes across maturity. Ceteris 
paribus, diminishing f3 means raising the barrier and thus diminishing volatilities, 
but a has a lower bound equal to zero. So, in the end, the feasible interval for the (3 
values is rather narrow and /? = 0.08 is one value in this narrow interval. For previous 
dates, the calibration was possible also with other values of /3 but we decided to keep 
it fixed to 0.08 for uniformity reasons. 
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Figure 7: Same as above on December 10th, 2003. 
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Figure 8: Volatility term structure T i— > a{T) on December 8th, 2003. 

From the results (Tables [TT] to [T4|) . we see that in this second framework, the 
stripped volatilities are much higher than before, especially in proximity of the cri- 
sis, where they are suspiciously high for practical purposes and raise a number of 
questions on realism. However, this effect is essentially due to the low values of H 
and not to (3, as one could suspect. Indeed, in Table [15] we present the results of 
a calibration performed for the first date, with f5 = 0.08 and H computed with the 
"credit spread based" method, and it is immediately clear that the volatilities are 
no longer pathologically large and are nearly the same as in Table [3 Anyway, the 
survival probabilities are not sensibly influenced by the method chosen to estimate 
H when the CDS to be calibrated are fixed; this lack of sensitivity is mainly due to 
the peculiar role of f3 and to the possibility to adjust its value, which in practice is 
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Figure 9: Survival probability on December 10th, 2003. 
an important degree of freedom of the model. 

Remark 4.2. (SBATIP, Scenario Barrier version of the ATIP model) A 

final remark on the possible use of random scenario-based barriers is in order. It 
is well known, see for example Bielecki and Rutkowski (2001), that in structural 
models with deterministic threshold barriers the default time is predictable. This is 
at times considered a drawback, especially in multi-name situations where one has to 
take into account contagion effects (Giesecke (2002)). What can help in this case is 
a random scenario-based default barrier. We may see that in our specific case the 
introduction of a random variable for the H parameter, independent of the driving 
Brownian Motion W , induces through iterated conditioning an analytical formula for 
default probabilities and thus CDS's. This formula is based on a one dimensional 
integration of (Q^j in H against the density of the random variable replacing the 
deterministic H. We call this model the SBATIP (Scenario Barrier Analytically 
Tractable 1st Passage model) and we investigate its use in Brigo and Tarenghi (2005), 
where we use parameters in the distribution of H to try and calibrate CDS prices. 
Scenarios on value of the firm volatilities can be employed as well (SVBATIP model). 

5 A fundamental example: Pricing Counterparty 
Risk in Equity Swaps 

In this section we present an example of pricing with the calibrated structural model. 
This example concerns the valuation of an equity swap where we take into account 
counterparty risk, and is chosen to highlight one case where the calibrated structural 
model may be preferable to a reduced-form intensity model calibrated to the same 
market information. This is an illustration of a more general situation that typically 
occurs when one tries to price the counterparty risk in an equity payoff. We will see 
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that it is possible to split the expectation of the payoff, and that the decomposition 
roughly involves the valuation of the same payoff without counterparty risk and the 
valuation of an option on the residual NPV of the considered payoff at the default time 
of the counterparty. Therefore including the counterparty risk adds an optionality 
level to the payoff. 

Let us consider an equity swap payoff. Assume we arc a company "A" entering 
a contract with company "B" , our counterparty. The reference underlying equity is 
company "C" . The contract, in its prototypical form, is built as follows. Companies 
"A" and "B" agree on a certain amount K of stocks of a reference entity "C" (with 
price S) to be taken as nominal (A^ = K Sq). The contract starts in = and 
has final maturity Ti, = T . At t = there is no exchange of cash (alternatively, we 
can think that "B" delivers to "A" an amount K of "C" stock and receives a cash 
amount equal to KSq). At intermediate times "A" pays to "B" the dividend flows of 
the stocks (if any) in exchange for a periodic rate (for example a semi-annual LIBOR 
or EURIBOR rate L) plus a spread X. At final maturity T = Tj,, "A" pays KSt to 
"B" (or gives back the amount K of stocks) and receives a payment KSq. This can 
be summarized as follows: 

Initial Time 0: no flows, or 
A — > KSo cash — > B 
A < — K equity < — B 

Time Tj: 
A — > equity dividends — > B 
A < — Libor + Spread < — B 

Final Time T;,: 
A — > K equity — > B 
A < — KSo cash < — B 

The price of this product can be derived using risk neutral valuation, and the 
(fair) spread is chosen in order to obtain a contract whose value at inception is zero. 
We ignore default of the underlying "C" , thus assuming it has a much stronger credit 
quahty than the counterparty "B" . This can be the case for example when "C" is an 
equity index (Pignatelli (2004)). It can be proved that if we do not consider default 
risk for "B" , the fair spread is identically equal to zero. But when taking into account 
counterparty default risk in the valuation the fair spread is no longer zero. In case 
an early default of the counterparty "B" occurs, the following happens. Let us call 
T — tb the default instant. Before r everything is as before, but if r < T, the net 
present value (NPV) of the position at time r is computed. If this NPV is negative 
for us, i.e. for "A", then its opposite is completely paid to "B" by us at time r itself. 
On the contrary, if it is positive for "A", it is not received completely but only a 
recovery fraction Rbc of that NPV is received by us. It is clear that to us ( "A" ) the 
counterparty risk is a problem when the NPV is large and positive, since in case "B" 
defaults we receive only a fraction of it. 

Analytically, the risk neutral expectation of the discounted payoff is {L{S, T) is 
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the simply compounded rate at time S for maturity T): 



r r ^ 

n^5(0) = Eo 1{.>T,} -i^NPV°f,2,„,,(0)+ir5o$^D(0,T,)a,(L(T,_i,T, 

L i=i 



+X)+D{0,n){KSo-KST,) 



+1 



/3(r)-l 



i=l 



+X) +D(0,r)(REc(NPV(r))+ - (-NPV(r))+) 



(22) 



where 



r ^ 

NPV(r) = K\-KNPV:^^^,^iT) + KSo D{r,T,)a,{L{T,_,,T,) + X) 

^ i=/3(r) 

+ {KSo-KSt,)D{t,T,)\. (23) 



We denote by NPV^^*^g^^^(M) the net present value of the dividend flows between s 
and t computed in u. 

In the following we will prove the 

Proposition 5.1. (Equity Return Swap price under Counterparty Risk). 

The fair price of the Equity Swap defined above, i.e. [2^) . can he simplified as follows: 

^ r 

IiEs{Q) = ir5oX5^a,P(0,T,) - Lao¥.A l{.<r,}P(0, r)(iVPy(r)) 
i=i 

The first term is the equity swap price in a default-free world, whereas the second one 
is the optional price component due to counterparty risk. 

Remark 5.2. (Different choice for the NPV). Let us focus for a second on the 
NPV term in [23^) . In particular let us concentrate on the interval [T^(t-)_i, T^(t-)) 
containing the default time r. If r is strictly inside the interval, the LIBOR + spread 
resetting at the earlier T^(^)_i is paid for the whole interval. 

A different formulation is the following. In case the default time is inside a certain 
interval, the LIBOR + spread term for the part of this interval preceding default is 
immediately paid and we add the remaining part, to he paid at the end of the interval, 
discounted hack from the end of the interval to the current (default) time. This 
amounts to let the summation in ^23) start from /9(r) + 1 rather than from [3{t) and 
to add the following term in the NPV payoff: 

{L {Tp(r)-i,Tp(T)) +X) (r-T^(^)_i) + (L [Tp(r)-i,Tp{T)) + X) [T^^^) - t) D {T,Tp(^r)) 

Formulation [2^) reads like postponing the default t to the first Ti following r and 
then discounting it hack to the current (default) time. 
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Both formulations make sense so contract specifications are important in order to 
price this kind of products. Analytical pricing formulas with our choice result to be 
more tractable than with the alternative choice presented in this remark, where small 
accruing terms appear. Such accruing terms are very small in general and can be 
ignored. 

If we try and find the above price by computing the expectation through a Monte 
Carlo simulation, we have to simulate both the behavior of the equity "C" underlying 
the swap, which wc call St = S^ , and the default of the counterparty "B". In 
particular we need to know exactly r = Tb- Obviously the correlation between "B" 
and "C" could have a relevant impact on the contract value. Here the structural model 
can be helpful: Suppose to cahbrate the underlying process V to CDS's for name "B" , 
finding the appropriate default barrier and volatilities according to the procedure 
outlined earlier in this paper with the ATIP model. We could set a correlation 
between the processes V^^ for "B" and St for "C", derived for example through 
historical estimation directly based on equity returns, and simulate the joint evolution 
of [^^^,5*^]. As a proxy of the correlation between these two quantities we may 
consider the correlation between Sf and S^ , i.e. between equities. 

Remark 5.3. (Using reduced form models?) At this point it may be interesting 
to consider the alternatives offered by an intensity reduced form model. If one takes 
a deterministic intensity, since Poisson processes and Brownian Motions defined on 

the same space are independent, there is no way to introduce a correlation between 
the default event for "B" and the value of the equity "C". To take into account this 
correlation we may have to move to stochastic intensity. However, correlating the 
stochastic intensity of "B" to the equity of "C" may lead to a poorly tractable model, 
even under tractable intensities as in Brigo and Alfonsi (2005), and the amount of 
correlation induced in this way may he too small for practical purposes. Furthermore, 
historical estimation of the correlation between the instantaneous intensity and the 
equity of "C" may pose some problems that the structural model avoids by definition. 

Now let us go more into details in analyzing the payoff. For simplicity, we as- 
sume deterministic dividends and rates, so that wc can substitute the discount factors 
D{t,T) with zero coupon bond prices P{t,T) and move them outside the expecta- 
tions. Moreover we recall that L{Ti_i,Ti) = (P(0, Ti_i)/P(0, T^) - 1) /a^ and, as 
a consequence, Ylf=A+iPi^^Ti)aiL{Ti_uTi) = P{0,Ta) - P{0,Tb). Finally, we set 
Rec = 1 — Lgd, consistently with CDS notation, and notice that 

Rbc (NPV(r))+ - (-NPV(t))+ = NPV(r) - Lgd(NPV(t))+. (24) 

We can rewrite the discounted NPV as 

b 

P(0,T)NPV(r) = -KNPYZjLndsiO) + KSo P{0,T,)a,{m.,,T,) + X) 
+KSoP{0,n) - KEr{STjP{0,n). (25) 
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Substituting and in ([22]) we obtain: 



UEsiO) = Eo<^l{.>T,} 



+P{0,n){KSo-KST,) 



+ l{r<T6} 



i=l 



/3(r)-l 



+P(0,r)NPV(r) -LGDP(0,r)(NPV(r))- 



i=l 



+P{Q,T,)KSo-P{Q,T,)KSt, 



+ l{r<T6} 



i=l 



+P{Q,T,)KSo - P{0,n)KE,{ST,} - LgdP(0, r)(NPV(r))- 



-i^NPV 



dividends ^ 



(0) + KSo J2 ^(0' Ti) + X) 

i=l 

+P(0,n)i^5o - P(0,Tfe)i^Eo|l{,>T,}^T,| - P(0,Tfe)irEo|l{.<T,}E4^Tj| 



LgdEo l{,<T,}P(0,r)(NPV(r))- 



-^ NPv°^,S,„,,(o) + i^^o + KSoX J2 t;, 



i=l 



P(0,T,)irEo i{r>n}Sn + i{r<n}Er{Sn} 



LGDEo<^l{.<T,}P(0,r)(NPV(r))- 



(26) 



Now, observing that 

Eo{l{r>T,}5'ri, + l{r<Ti,}IEr{5'Tj} 



Eo{l{r>n}Sn + (1 - 1{.>t,})E,{StJ} 
Eo{E.{5rJ + l{r>n}{Sn - K{Sn})} 
Eo{e.{5tJ} + MUr>n}iSn - Sn)} 
So 



P(0,T,) P(0,Tfc) 



(27) 
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the expected payoff (12^ becomes 

b 

i=l 

^m^J l^P(0,T,) P(0,T, 
-LGDEo|l{.<T}P(0,r)(NPV(r)) 



= ir^oXj]a,P(0,T,)-LGDEo|l{.<T,}P(0,r)(NPV(r))+|. (28) 

which proves our proposition above. As noticed in the proposition, the first term 
in the last expression in (|28l) is n|;5r^^^(0), i.e. the payoff of the Equity Swap in a 
default free ("defree") world. Indeed, as stated before, in a default free market the 
fair spread would be zero, so that if we add a nonzero spread X the value of the 
equity swap in a default free world is simply the NPV of the future spread cash flows. 

The second term is very important if we allow counterparty "B" to default. As it 
is non-positive, to have a fair value for n£;5(0) we need a non-negative spread X. 

At this point it may be worth noticing that the above decomposition is rather 
general. When we include counterparty risk into the valuation, we always obtain a 
decomposition 

Price = Priceo cfaultFrec - (1-recovery) CallOption on residual NPV 

The call option has strike 0. As we have observed earlier, the counterparty risk 
introduces, among other features, a further optionality level into the valuation. 

Going back to our equity swap, now it is possible to run the Monte Carlo simu- 
lation, looking for the spread X that makes the contract fair. The simulation itself 
is simpler when taking into account the following computation included in the dis- 
counted NPV: 

P{T,n)E^{ST,} =Sr- ^Py^Sendsir) (29) 

SO that we have 

b 

P(0,r)NPV(r) = KSo P(0, T,)«,,(L(T,_i, T,) + X) + iT^om ^fe) 

i=l3{T) 

-KP{0,T)Sr = 
b 

= KSo PiO,T,)a,X + KSoPiO,Tpir)-i) 

i=l3{T) 

-KP{0,T)Sr. (30) 

The reformulation of the original expected payoff fl22l) as in fl28l) presents an 
important advantage in terms of numerical simulation. In fact in fl2^ we have a 
global expectation, hence we have to simulate the exact payoff for each path. In fl25]) . 
with many simplifications, we have isolated a part of the expected payoff out of the 
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main expectation. This isolated part has an expected value that we have been able to 
calculate, so that it does not have to be simulated. Simulating only the residual part 
is helpful because now the variance of the part of the payoff that has been computed 
analytically is no longer affecting the standard error of our Monte Carlo simulation. 
The standard error is indeed much lower when simulating fl25]) instead of ([22D- The 
expected value we computed analytically above involves terms in St^ which would 
add a lot of variance to the final payoff. In (!28|) the only St^ term left is in the 
optional NPV part. 

We performed some simulations under different assumptions on the correlation 
between "B" and "C". We considered five cases: p = —1, p = —0.2, p = 0, p = 0.5 
and p = 1. In Tabled] we present the results of the simulation, together with the error 
given by one standard deviation (Monte Carlo standard error). For counterparty "B" 
we used the same CDS rates of the company analyzed in Section [31 For the reference 
stock "C" we used a hypothetical stock with initial price 5*0 = 20, volatility a = 20% 
and constant dividend yield q = 0.80%. The contract has maturity T = 5y and the 
settlement of the LIBOR rate has a semi-annual frequency. Finally, we included a 
recovery rate Rec = 40%. The starting date is the same we used for the calibration, 
i.e. March 10th, 2004. Since the reference number of stocks K is just a constant 
multiplying the whole payoff, without losing generality we set it equal to one. 

In order to reduce the errors of the simulations, we have adopted a variance 
reduction technique using the default indicator (whose expected value is the known 
default probability) as a control variate. In particular we have used the default 
indicator l{r<r} at the maturity T of the contract, which has a large correlation with 
the final payoff. Even so, a large number of scenarios are needed to obtain errors with 
a lower order of magnitude than X. In our simulations we have used = 2000000. 

We notice that X increases together with p. This fact can be explained in the 
following way. Let us consider the case of positive correlation between "B" and "C" : 
This means that, in general, if the firm value for "B" increases, moving away from 
the default barrier, also the stock price for "C" tends to increase due to the positive 
correlation. Both processes will then have high values. Instead, again under positive 
correlation, if V^^ lowers towards the default barrier, also will tend to do so, going 
possibly below the initial value 5*0. In this case NPV{t) has a large probability to be 
positive (see fl30l) ). so that one needs a large X to balance it, as is clear when looking 
at the final payoff fl28|) . On the contrary, for negative correlation, the same reasoning 
can be applied, but now if V^^ lowers and tends to the default barrier, in general 

will tend to move in the opposite direction and the corresponding NPV{t) will 
probably be negative, or, if positive, not very large. Hence the "balancing" spread 
X we need will be quite small. 

The simulation error depends on p too, and in particular it is influenced by two 
main effects. The first one is the same as before. Indeed, a positive correlation means 
high probability of positive NPV in case of default. As the default probability of "B" 
does not depend on p, it is clear that a positive correlation leads to evaluate the 
expectation of a random variable which has a high probability of being different from 
zero (and in particular positive), while a negative correlation leads to the expectation 
of a variable which is almost always zero, leading to a lower Monte Carlo error. The 
other effect is due to the impact of p on the correlation between the payoff variate 
and the control variate. It happens that this correlation is decreasing in p: It is 
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p 


A 


tib payoii 


MC error 


1 

-i 


n 
u 


n 
u 


n* 
u 


-0.2 


2.45 


-0.02 


1.71 





4.87 


-0.90 


2.32 


0.5 


14.2 


-0.53 


2.71 


1 


24.4 


-0.34 


0.72 



Table 6: Fair spread X (in basis points) of the Equity Swap in five different correlation 
cases, 5*0 = 20. We also report the value of the average of the simulated payoff (times 10000) 
across the 2000000 scenarios and its standard error, for the chosen X, thus showing that 
X is indeed fair in practice, since all the averages and MC windows are close to zero. X's 
have been found by iterating the MC simulation for different values of X until the payoff 
was found to be sufficiently small. The MC error for p = —1 (marked with an asterisk) is 
null since in each simulated scenario the simulated term has been found to be identically 
zero (i.e. the NPV was always negative). 

nearly —1 when p = 1, and then increases with lower values of p (it is about —0.4 for 
p = -0.2). 

As a final remark we notice that the spread X corresponding to p = — 1 is iden- 
tically equal to zero. The reason is quite simple: When we have maximum negative 
correlation between the counterparty "B" and the reference entity "C", the NPV in 
( l28l) is negative in all the simulated scenarios involving default, so that the related 
term in the price of the equity swap is always zero. We are then left only with the 
fixed term and this is null only when X is zero. All this reasoning also explains why 
the error of the simulation, marked with an asterisk, is zero. 

6 Conclusions 

In general the link between default probabilities and credit spreads is best described 
by intensity models. The credit spread to be added to the risk free rate represents 
a good measure of a bond credit risk for example. Yet, intensity models present 
some drawbacks: They do not link the default event to the economy but rather to 
an exogenous jump process whose jump component remains unexplained from an 
economic point of view. Moreover, when dealing with default correlation, a copula 
function must be introduced between the jump processes thresholds in a way that 
has no clear immediate relation with equity correlation, the only source of correlation 
that can be used for practical purposes. 

In this paper we introduced an analytically tractable structural model that allows 
for a solution to the above points. In this model the default has an economic cause, 
in that it is caused by the value of the firm hitting the safety barrier value, and 
all quantities are basic market observables. Also, when dealing with multi-name 
products, the model allows for the introduction of the correlation in a very natural 
way, by simply correlating shocks in the different values of the firm equities by means 
of equity correlation. 

We showed how to calibrate the model parameters to actual market data: Start- 
ing from CDS quotes, we calibrated the value of the firm volatilities and found the 
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barrier triggering the default that was consistent with CDS quotes and also leading 
to analytical tractability. We also explained the analogies with barrier option pricing, 
in particular the case with time dependent parameters. We hinted at the possible 
use of a scenario based random barrier, as in Brigo and Tarenghi (2005). 

As a practical example, we also applied the model to a concrete case, showing 
how it can describe the proximity of default when time changes and the market 
quotes for the CDS's reflect increasing credit deterioration. When the market detects 
a company crisis, it responds with high CDS quotes and this translates into high 
default probabilities, i.e. high probabilities for the underlying process to hit the 
safety barrier, that in turn translate in high calibrated volatilities for the firm value 
dynamics. We also took advantage of the numerical examples to further explain and 
interpret a few model parameters whose role needed further investigation. 

The calibrated model is useful for pricing more sophisticated derivatives depend- 
ing on default. We gave the fundamental example of equity swaps pricing under 
counterparty risk. 

7 Appendix: Results of the case study calibrations 
and market data 

7.1 Credit spread based barrier parameter H 

Here we present the numerical results of the calibration relative to the case study 
discussed in Section HI For each date we present the spot price of the equity with its 
historical volatility, the rates and the recovery rates of the (running) CDS's (quarterly 
paid), the outputs of the calibration. For all cases we used P = 0.5. H is computed 
using the credit spread based method explained in Section [3l 

• September 10th, 2003. 

Stock Price = 2.898 euros, volatility = 5%, recovery rate = 40%. The barrier 
parameter is H = 0.8977 (re-scaled Vq = 1). 



CDS maturity 


Rate (bps) 


Volatility nodes 


Survival 


lO-sep-03 




5.012% 


100.000% 


ly 


192.5 


5.012% 


96.673% 


3y 


215 


3.103% 


89.524% 


5y 


225 


3.178% 


82.471% 


7y 


235 


3.551% 


75.375% 


lOy 


235 


3.658% 


66.998% 



Table 7: Outputs of the calibration on September 10th, 2003. Here /3 = 0.5 and H = 
0.8977. 

• November 28th, 2003. 

Stock Price = 2.297 euros, volatility = 14%, recovery rate = 40%. The barrier 
parameter is H = 0.8052 (re-scaled Vq = 1). 
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CDS maturity 


Rate (bps) 


Volatility nodes 


ouivival 


28-nov-Ud 




"1/1 o 1 (7/ 

14.08170 


1 AA AAAQ/ 

lOO.OOO/o 




1 zo 


l^.Uol /O 


O / .DZU/0 


3y 


630 


10.800% 


72.529% 


5y 


570 


11.489% 


62.561% 


7y 


570 


16.235% 


51.810% 


lOy 


570 


22.793% 


39.217% 



Table 8: Outputs of the calibration on November 28th, 2003. Here /3 = 0.5 and H = 0.8052. 

• December 8th, 2003. 

Stock Price = 2.237 euros, volatility = 20%, recovery rate = 25%. The barrier 
parameter is if = 0.7730 (re-scaled Vq — I). 



CDS maturity 


Rate (bps) 


Volatility nodes 


Survival 


8-dec-03 




20.197% 


100.000% 


ly 


1450 


20.197% 


81.289% 


3y 


1200 


17.972% 


62.065% 


5y 


940 


13.685% 


56.263% 


7y 


850 


18.771% 


49.386% 


lOy 


850 


36.661% 


35.626% 



Table 9: Outputs of the calibration on December 8th, 2003. Here /3 = 0.5 and H = 0.7730. 

• December 10th, 2003. 

Stock Price = 2.237 euros, volatility = 50%, recovery rate — 15%. The barrier 
parameter is H = 0.7253 (re-scaled Vq = 1). 



CDS maturity 


Rate (bps) 


Volatility nodes 


Survival 


lO-dec-03 




50.000% 


100.000% 


ly 


5050 


50.000% 


55.452% 


3y 


2100 


4.325% 


55.208% 


5y 


1500 


19.950% 


50.910% 


7y 


1250 


24.063% 


46.705% 


lOy 


1100 


37.422% 


39.121% 



Table 10: Outputs of the calibration on December 10th, 2003. Here /? = 0.5 and H = 
0.7253. 

Remcirk 7.1. In this last case the calibration algorithm for the deterministic 
intensity model could not achieve a solution, so it was not possible to extract 
the suitable value for the "credit spread based" H. In such cases we can proceed 
in an alternative way. Let a he the equity volatility for the first maturity that 
is given by the equity market. 



a. Choose a first guess H^^^ for H. Set i—1; 



D. Brigo, M. Tarenghi: CDS Calibrated Structural Model and Application to Equity Swaps 32 



b. Calibrate the structural model to the first maturity CDS market quote 
through the parameter a when the barrier is set to H^'^\- Let cr*^*^ be the 
calibrated volatility. 

c. Compute the first maturity default probability with the structural model 
characterized by H^'^^ and cr^*^ . 

0?. Solve in H the equation where we equate the default probability computed 
in "c" above to the default probability of a structural model characterized 
by H and a; Set to the found value of H; Replace i by i + 1 and 

restart from "b". 

Typically, after a number of iterations, o"^'^ will be very close to a and we may 
set H to the corresponding value if . 

7.2 Excursion/protection analogy barrier parameter H 

Here we present the results of the cahbration performed using the excursion protection 
analogy to estimate H, i.e. H/Vq — Rec. Now (3 — 0.08 and the other parameters 
are the same as before. 

• September 10th, 2003. 

The barrier parameter is if = Rec = 0.4 (re-scaled Vq = 1). 



CDS maturity 


Rate (bps) 


Volatihty nodes 


Survival 


lO-sep-03 




42.713% 


100.000% 


ly 


192.5 


42.713% 


96.674% 


3y 


215 


26.588% 


89.524% 


5y 


225 


27.343% 


82.473% 


7y 


235 


30.686% 


75.377% 


lOy 


235 


31.778% 


67.002% 



Table 11: Outputs of the calibration on September 10th, 2003. Here /3 = 0.08 and H = 0.4. 

• November 28th, 2003. 

The barrier parameter is H — Rec = 0.4 (re-scaled Vq — I). 



CDS maturity 


Rate (bps) 


Volatility nodes 


Survival 


28-nov-03 




58.808% 


100.000% 


ly 


725 


58.808% 


87.613% 


3y 


630 


44.225% 


72.512% 


5y 


570 


46.237% 


62.541% 


7y 


570 


63.898% 


51.780% 


lOy 


570 


85.758% 


39.146% 



Table 12: Outputs of the calibration on November 28th, 2003. Here /3 = 0.08 and H = 0.4. 
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• December 8th, 2003. 

The barrier parameter is H = Rec = 0.25 (re-scaled Vq = 1). 



CDS maturity 


Rate (bps) 


Volatihty nodes 


Survival 


8-dec-03 




107.846% 


100.000% 


ly 


1450 


107.846% 


81.277% 


3y 


1200 


94.506% 


62.041% 


5y 


940 


71.046% 


56.246% 


7y 


850 


96.431% 


49.369% 


lOy 


850 


181.575% 


35.559% 



Table 13: Outputs of the cahbration on December 8th, 2003. Here /3 = 0.08 and H = 0.25. 

• December 10th, 2003. 

The barrier parameter is H = Rec = 0.15 (re-scaled Vq = 1). 



CDS maturity 


Rate (bps) 


Volatility nodes 


Survival 


lO-dec-03 




292.060% 


100.000% 


ly 


5050 


292.060% 


55.381% 


3y 


2100 


21.603% 


55.198% 


5y 


1500 


114.713% 


50.904% 


7y 


1250 


137.488% 


46.700% 


lOy 


1100 


210.474% 


39.104% 



Table 14: Outputs of the cahbration on December 10th, 2003. Here /3 = 0.08 and H = 0.15. 

Finally, we present the calibration performed for the first date (September 10th, 2003) 
with l3 = 0.08 and H estimated by means of the credit spread method. Comparing 
the values in Table [TS] with those in Table [7] we see that the volatilities are nearly the 
same, and this holds in general, at least when the volatilities are not too high. Here 
the barrier parameter is H = 0.8969 (re-scaled Vq = 1). 



CDS maturity 


Rate (bps) 


Volatility nodes 


Survival 


lO-sep-03 




5.012% 


100.000% 


ly 


192.5 


5.012% 


96.673% 


3y 


215 


3.064% 


89.522% 


5y 


225 


3.110% 


82.468% 


7y 


235 


3.445% 


75.369% 


lOy 


235 


3.508% 


66.988% 



Table 15: Outputs of the cahbration on September 10th, 2003. Here (3 = 0.08 and H = 
0.8969. 
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7.3 Discount Curves 

In the following we report the Discount Curves that have been used in the simulations. 



• September 10th, 2003. 



Date 


Discount 


Date 


Discount 


Date 


Discount 


Sep 10th, 03 


1.00000 


Jun 14th, 04 


0.98320 


Sep 12th, 


11 


0.71837 


Sep 11th, 03 


0.99994 


Sep 13th, 04 


0.97721 


Sep 12th, 


12 


0.68171 


Sep 15th, 03 


0.99971 


Sep 12th, 05 


0.94740 


Sep 12th, 


13 


0.64624 


Sep 19th, 03 


0.99948 


Sep 12th, 06 


0.91195 


Sep 12th, 


14 


0.61204 


Oct 13th, 03 


0.99805 


Sep 12th, 07 


0.87373 


Sep 14th, 


15 


0.58032 


Nov 12th, 03 


0.99627 


Sep 12th, 08 


0.83501 


Sep 12th, 


16 


0.54942 


Dec 12th, 03 


0.99447 


Sep 14th, 09 


0.79557 


Sep 12th, 


17 


0.51882 


Mar 12th, 04 


0.98896 


Sep 13th, 10 


0.75666 


Sep 12th, 


18 


0.48937 



Table 16: Discount curve on September 10th, 2003. 



• November 28th, 2003. 



Date 


Discount 


Date 


Discount 


Date 


Discount 


Nov 28th, 03 


1.00000 


Sep 02nd, 04 


0.98200 


Dec 02nd, 11 


0.70648 


Dec 01st, 03 


0.99983 


Dec 02nd, 04 


0.97508 


Dec 03rd, 12 


0.66950 


Dec 03rd, 03 


0.99971 


Dec 02nd, 05 


0.94202 


Dec 02nd, 13 


0.63449 


Dec 09th, 03 


0.99936 


Dec 04th, 06 


0.90424 


Dec 02nd, 14 


0.60081 


Jan 02nd, 04 


0.99792 


Dec 03rd, 07 


0.86449 


Dec 02nd, 15 


0.56852 


Feb 02nd, 04 


0.99607 


Dec 02nd, 08 


0.82431 


Dec 02nd, 16 


0.53760 


Mar 02nd, 04 


0.99435 


Dec 02nd, 09 


0.78427 


Dec 04th, 17 


0.50786 


Jun 02nd, 04 


0.98849 


Dec 02nd, 10 


0.74490 


Dec 03rd, 18 


0.47965 



Table 17: Discount curve on November 28th, 2003. 



• December 8th, 2003. 



Date 


Discount 


Date 


Discount 


Date 


Discount 


Dec 08th, 03 


1.00000 


Sep 10th, 04 


0.98271 


Dec 12th, 


11 


0.71134 


Dec 09th, 03 


0.99994 


Dec 10th, 04 


0.97612 


Dec 10th, 


12 


0.67460 


Dec nth, 03 


0.99983 


Dec 12th, 05 


0.94456 


Dec 10th, 


13 


0.63960 


Dec 17th, 03 


0.99949 


Dec nth, 06 


0.90787 


Dec 10th, 


14 


0.60587 


Jan 12th, 04 


0.99793 


Dec 10th, 07 


0.86880 


Dec 10th, 


15 


0.57341 


Feb 10th, 04 


0.99620 


Dec 10th, 08 


0.82911 


Dec 12th, 


16 


0.54209 


Mar 10th, 04 


0.99448 


Dec 10th, 09 


0.78937 


Dec nth. 


17 


0.51229 


Jun 10th, 04 


0.98882 


Dec 10th, 10 


0.74996 


Dec 10th, 


18 


0.48373 



Table 18: Discount curve on December 8th, 2003. 
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• December 10th, 2003. 



Date 


1 )m( -ni lilt 


Date 


1 )i <-;(■( unit 


Date 


Flmcrnmt 

i_y iov. y J liiiu 


Dpp 1 nth n*^ 


1 nnnnn 


Sen 1 '^th 04 




Dec 12th 


11 


71 1 65 


Dec nth, 03 


0.99994 


Dec 13th, 04 


0.97610 


Dec 12th, 


12 


0.67476 


Dec 15th, 03 


0.99972 


Dec 12th, 05 


0.94476 


Dec 12th, 


13 


0.63957 


Dec 19th, 03 


0.99949 


Dec 12th, 06 


0.90802 


Dec 12th, 


14 


0.60571 


Jan 12th, 04 


0.99805 


Dec 12th, 07 


0.86889 


Dec 14th, 


15 


0.57303 


Feb 12th, 04 


0.99621 


Dec 12th, 08 


0.82922 


Dec 12th, 


16 


0.54191 


Mar 12th, 04 


0.99448 


Dec 14th, 09 


0.78919 


Dec 12th, 


17 


0.51196 


Jun 14th, 04 


0.98874 


Dec 13th, 10 


0.74987 


Dec 12th, 


18 


0.48324 




Table 19: Discount curve on December 10th, 2003. 






• March 10th, 2004. 












Date 


Discount 


Date 


Discount 


Date 


T^isrnnnt, 


Mar 10th, 04 


1.00000 


Dec 13th, 04 


0.98422 


Mar 12th, 


12 




Mar nth, 04 


0.99994 


Mar 14th, 05 


0.97884 


Mar 12th, 


13 


0.69831 


Mar 15th, 04 


0.99972 


Mar 13th, 06 


0.95375 


Mar 12th, 


14 


0.66271 


Mar 19th, 04 


0.99949 


Mar 12th, 07 


0.92279 


Mar 12th, 


15 


0.62817 


Apr 13th, 04 


0.99807 


Mar 12th, 08 


0.88770 


Mar 14th, 


16 


0.59552 


May 12th, 04 


0.99642 


Mar 12th, 09 


0.85063 


Mar 13th, 


17 


0.56392 


Jun 14th, 04 


0.99455 


Mar 12th, 10 


0.81234 


Mar 12th, 


18 


0.53344 


Sep 13th, 04 


0.98940 


Mar 14th, 11 


0.77355 


Mar 12th, 


19 


0.50440 



Table 20: Discount curve on March 10th, 2004. 
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